Spectral properties of quantized barrier billiards 
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The properties of energy levels in a family of classically pseudointegrable systems, the barrier 
billiards, are investigated. An extensive numerical study of nearest-neighbor spacing distributions, 
next-to-nearest spacing distributions, number variances, spectral form factors, and the level dy- 
namics is carried out. For a special member of the billiard family, the form factor is calculated 
analytically for small arguments in the diagonal approximation. All results together are consistent 
with the so-called semi-Poisson statistics. 
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I. INTRODUCTION 

Two decades after the first investigation of the quan- 
tum mechanics of nonintegrable polygonal billiards [0 , a 
renewed interest in this peculiar class of dynamical sys- 
tems has shown up recently. One reason is the fabrica- 
tion of polygonal-shaped optical microcavities |3[ ^, || . 
Another reason is the finding ^, |^ that certain planar ra- 
tional polygons (all angles between sides are of the form 
pTr/q, where p, q are relatively prime integers) have spec- 
tral properties very similar to those of mesoscopic disor- 
dered systems at the critical point of the metal-insulator 
transition [p| and to those of systems with interacting 
electrons |^ 

The classical dynamics in rational polygons having at 
least one critical corner with p > 1 is characterized as 
pseudointegrable . The phase space is foliated by two- 
dimensional invariant surfaces |l^ , like in integrable 
systems , but the genus of the surfaces is larger than 
one The flow on these surfaces is typically ergodic 
and not mixing JTsf . 

Pseudointegrable systems cannot be quantized accord- 
ing to the semiclassical Einstein-Brillouin-Keller rule [Q . 
As a consequence, the statistical properties of energy 
levels of classically pseudointegrable systems are differ- 
ent from those of integrable systems which are generi- 
cally well described by Poissonian random processes . 
For example, the nearest-neighbor spacing distribution 
of pseudointegrable systems generically displays a clear 
level repulsion |^ as in the case of the Gaussian orthogo- 
nal ensemble (GOE) of random-matrix theory which 
describes fully chaotic systems with time-reversal sym- 
metry Significant deviations from GOE are ob- 
served first theoretically ]l7t and later experimentally in 
microwave cavities . 

Recently, the semi-Poisson (SP) statistics have been 
proposed as reference point for the spectral statistics of 
pseudointegrable systems Following ||^, |l^, we 

define SP statistics by removing every other level from 
an ordered Poisson sequence {a;„}. Unfortunately, the 
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term "semi-Poisson" was originally coined for a sequence 
{Un} where ?/„ = (x„ -I- a;„+i)/2 |6|. We call this here 
interpolated-Poisson (IP) statistics. IP and SP statistics 
have identical nearest-neighbor spacing distributions but 
other spectral quantities in general differ. 

The SP conjecture has been verified numerically for 
right triangular |^ and rhombus billiards where 
only small differences to SP have been found. However, 
numerical works on right triangles in a regime of high 
level numbers seem to indicate that the statistical prop- 
erties are nonstationary |2^, ^ with increasing, but still 
small, deviations from SP as the energy is increased . 
For certain right triangles, it has been shown analyti- 
cally that the spectral form factor for small arguments is 
located around the corresponding SP result p3[ . 

However, the triangles studied in are not generic 
rational polygons, because they belong to the class of 
Veech polygons which may have special spectral 
properties. In this paper, we study the symmetric bar- 
rier billiards ||2|, ||, |^, ||, H, |3^ where the even- 
symmetry states, the "pure barrier-billiard states" , are 
expected to show the generic behavior. We provide ana- 
lytical and extensive numerical calculations showing that 
the spectral properties are fully consistent with SP statis- 
tics. Moreover, our results throw some light on the non- 
stationarity observed in ||2l|, . 

The paper is organized as follows. After defining the 
billiard family in Sec. II, we compute analytically the 
form factor for small arguments in Sec. pil|. Numerical 
results are presented in Sec. [V, Section M contains a 
conclusion. 



II. BARRIER BILLIARDS 

The family of barrier billiards consists of rectangles 
with sizes Z^, ly and a barrier placed on the symmetry line 
X = Ix/^ as shown in Fig. |l|(a). The length of the barrier 
I S (0, ly) is the only nontrivial parameter. The free 
motion of a point particle with mass m and momentum 
iPx,Py) bounded by elastic refiections at the boundary of 
the billiard has a second constant of motion K — p^ in 
addition to Hamilton's function H. Hence, the dynamics 
in phase space is restricted to invariant surfaces {H, K) = 
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constant. The topology of these surfaces is not that of 
a torus (with genus 1) but that of a two-handled sphere 
(genus 2) due to the critical corner at the upper end of the 
barrier; see for the relation between critical corners 
and the genus of invariant surfaces in pseudointegrable 
billiards. 

The billiards are Veech (roughly speaking, this prop- 
erty implies a special kind of hidden symmetry) if and 
only if Z/Zy is a rational number Still, a typical sym- 
metric barrier billiard is not a generic pseudointegrable 
system since it is composed of two copies of an integrable 
sub-billiard, the rectangle shown in Fig. |l|(b). This prop- 
erty is identical to almost-integrability |32| in the case of 
l/ly being rational. 
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FIG. 1: (a) Trajectory (dotted) in the full barrier billiard, a rectan- 
gle with a barrier between the points (x, y) = {Ix/'i, 0) and (ia;/2, 1). 
Symmetry-reduced system with (b) Dirichlet boundary conditions 
and (c) mixed boundary conditions: Dirichlet (Neumann) on solid 
(dashed) lines. 

The energy eigenstates are solutions of the Helmholtz 
equation with Dirichlet boundary conditions, i.e. van- 
ishing amplitude, on the boundary of the polygon. The 
states are odd or even with respect to the symmetry line. 
The former ones are trivial eigenstates of the integrable 
sub-billiard in Fig. [|(b). We therefore deal mainly with 
the even ones, the "pure barrier-billiard states", which 
fulfill mixed boundary conditions on the boundary of 
the symmetry-reduced polygon; see Fig. |l|(c). We expect 
that the pure barrier-billiard states show the generic fea- 
tures of energy states in rational polygons. 



III. THE SPECTRAL FORM FACTOR 

We here compute analytically a spectral quantity, the 
2-point correlation form factor, for the energy levels {En\ 
of a special member of the barrier-billiard family. Our 
calculation is inspired by that for the triangular billiards 
in p^ . It turns out that the present calculation is much 
simpler. As in ||2^ , we will apply the modern semiclas- 
sical theory based on trace formulas which express the 
density of states of a quantum system in terms of pe- 
riodic orbits of the underlying classical system |Q . For 
billiards, the semiclassical limit 7i ^ corresponds to the 
high-energy limit E ^ oo. Throughout the paper we use 
natural units such that T% — 2m — 1. 

The density of states can be written as sum of a smooth 



part and an oscillatory part 

d{E) =J2SiE- En) = d{E) + dosc{E). (1) 

n 

The fluctuations in the oscillatory part can be studied 
with the help of the 2-point correlation function 



R{€) = (dosc{E + ^)dosc{E ~ ^)) 



(2) 



Brackets (. . .) denote an energy averaging around E on 
an energy window much larger than the mean level spac- 
ing 1 /d{E), and much smaller than E. The Fourier trans- 
form of R{e) is the spectral form factor 
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We will concentrate on the limit 
Poisson i 



0; 



(3) 



1 for 



K(0) = 

_ 1/2 for SP l23| and for GOE @ ||. 
For two-dimensional rational polygons, the smooth 
part of the density of states is semiclassically described 
by Weyl's law d = A /{Air) where A is the area of the 
polygon; the oscillating part splits into two parts [p5| |3Q| 



dosciE) = dp.oXE) + d^^^XE) 



(4) 



The periodic orbit contribution 



dp o(E)=Y ^_^=e*^'''-*'^''"/2-»^/4 ^ ^^^^ 
^47ry2^ 

is a summation over classical (primitive and non- 
primitive) periodic orbits. These orbits are marginally 
stable and appear always in one-parameter families re- 
flecting the foliation of phase space by two-dimensional 
invariant surfaces, p labels these families; Ap denotes the 
surface in configuration space covered by a given family 
(without repetitions of primitive periodic orbits); Ip is 
the (non-primitive) length of periodic orbits; the Maslov 
index Vp is here twice the number of reflections at Dirich- 
let boundaries (Neumann boundaries do not contribute); 
k — y/E is the wave number. 

The diffractive orbit contribution dj o. (^) is a summa- 
tion over orbits starting and ending at critical corners of 
the polygon. This summation is more involved than the 
periodic orbit contribution ^T\ . In the limit r — > 0, how- 
ever, the form factor K{t) docs not depend on diffrac- 
tive orbits ||2^. With this insight a formula for K{0) 
has been derived in p3[ by inserting the periodic orbit 
contribution (|^) into Eq. (^ and employing the diagonal 
approximation (which is expected to be valid for small 
r) yielding 



K{0) = lim 



1 



SiT^d 



E 

P.O. 



I A 



-glSilp - Ankdr) , (6) 



where gp is the multiplicity of a given periodic-orbit fam- 
ily, i.e. the number of families with exactly the same 
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lengths, and the summation is performed over families 
with different lengths. 

For later considerations it is helpful to repeat the eval- 
uation of Eq. for the simplest case, the rectangular 
billiard, as done in |^^. A family of periodic orbits in a 
rectangle with sizes a, b {a and b are irrationally related) 
and area A = ab can be specified by two non-negative 
integers rup and rip, denoting the number of traversals 
across the billiard in the x- and y-direction, respectively. 
The length of each orbit is 



Ip = J{2mpaY + {2npbf 



(7) 



The number of periodic orbits N{1) up to length / is the 
number of lattice points in the positive (mp, np)-quadrant 
inside the ellipse (j^) asymptotically giving by 



N{1) 



16A 



(8) 



Due to the fact that all families cover the same area Ap = 
2A and have typically the same multiplicity Qp = 2 (time- 
reversal symmetry) the sum (|^) can be replaced by the 
following simple integral 
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if(0)= hm— = 



-S(l - 4TrkdT)^^:^^dl (9) 
I dl 



which gives -ft'(O) — 1 as expected for generic integrable 
systems Q. 

We now extend the previous calculation to the barrier 
billiard. To keep the calculation elementary, we restrict 
ourself to the special Veech case I = ly/2. The odd states 
are eigenstates in the rectangle with width a = lx/2^ 
height b = ly, and with Dirichlet boundary conditions. 



see Fig. |l|(b), so we get K{Q) = 1 as demonstrated above. 
The even states, the "pure barrier-billiard states" , fulfill 
mixed boundary conditions as shown in Fig. |l|(c). In 
the semiclassical trace formula (H), the inhomogeneous 
boundary conditions only influence the Maslov indices 
of the periodic orbits: a reflection at a Dirichlet bound- 
ary increases the index by two in contrast to a reflec- 
tion at a Neumann boundary which does not change the 
index. The resulting phase difference of tt between tra- 
jectories has an analog in billiards with a magnetic flux 
line where trajectories encircling a flux of 1/2 (in 
natural units) once pickup a phase tt. 

First, let us consider periodic orbits with fixed 
mp, Up > and nip odd. We write rUp = mN + rriD where 
niN^mD > count the number of reflections at x = lx/2 
with Neumann or Dirichlet boundary condition, respec- 
tively. Two cases have to be distinguished: m jv even and 
vtld odd; mjv odd and m/j even. The corresponding two 
types of orbits are related by a symmetry transformation 
(ignoring the boundary conditions), the reflection at the 
line y — ly/2. Hence, both types have the same Zp, gp 
and Ap. However, the Maslov indices are different due to 
the inhomogeneous boundary conditions: Vp mod 4 = 
and UpTooAA = 2, respectively. This implies that the 



contribution of both families to the trace formula (g) are 
identical differing just by a sign. Therefore both contri- 
butions cancel each other. 

Second, let us turn to periodic orbits with mp even. 
We begin with unfolding the orbits into a larger rect- 
angle with width a' = 2a. Assume, for simplicity, 
that nip = nip/2 is odd. The case m^ even can be 
treated by further unfolding of the orbits. Again, for 
ni'p — ni'j^ + m'jy odd there exists two kinds of periodic 
orbits related by symmetry: one with m'^ even and m'^, 



odd and one with 
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odd and 



even. These two 



kinds of trajectories either become congruent or remain 
separated when folded back into the original rectangle 
with width a. In the first case, we have to add the two 
different values of m^ and the two values of m'^ lead- 
ing to m^r and mu odd throughout the family. In the 
other case, we get m^r = 2m^ and m/^ = 2m'ij even 
since each orbit is symmetric with respect to the folding 
axis. Clearly, for fixed np only one of these two cases is 
possible. Hence, no cancellation occurs for even mp in 
the trace formula (|^) , in contrast to the complete cancel- 
lation in the case of nip odd. The simple consequence of 
which is that the number of periodic-orbit families which 
contributes to the trace formula (|^) is reduced by a fac- 
tor two. The same is true for the sum (|). From Eq. (|) 
follows then directly our main analytical result 



KiO) = - 



(10) 



Our calculated K{0) is not only close to the SP prediction 
as in the case of Veech triangles it agrees exactly 
with the SP prediction. 

The calculation for general barrier length is consid- 
erably more complicated. Yet, it should be possible to 
compute K{0) also for rational l/ly ^ 1/2 using methods 
developed in H 



IV. NUMERICAL RESULTS 

We here present numerical results on several statis- 
tical quantities for general symmetric barrier billiards. 
As representatives we choose the Veech billiard with 
I = ly/2 and one which is not Veech with I — lyUj, 
where ui = {Vb — l)/2 is (the reciprocal of) the golden 
mean. Irrationally related parameters Ix = 7r-\/87r/3 and 
ly = S-v/Stt/tt are taken. Billiards with I « and I « ly 
are not investigated since the semiclassical behavior of 
these limiting cases is expected to set in at extremely 
high energies. We consider two different energy regimes: 
(i) the medium-energy regime starting with the 40 000th 
level and ending with the 60 000th level, and (ii) the high- 
energy regime starting with the 400 000th level and end- 
ing with the 420 000th level. Our high-energy regime is 
below that of Ref. H and above that of Ref. ||2l|. 

We compute the eigenvalues with the mode-matching 
technique which is very efficient for barrier billiards as 
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described in detail in An accuracy of about 10~^ of 
the mean level spacing is achieved. 

To distinguish between local fluctuations in the level 
sequence Ei < E2 < E3 < . . . and a systematic global 
energy dependence of the average density we "unfold" 
the spectra in the usual way by setting En = N(En); 
see, e.g., N{E) is the smooth part of the integrated 
density of states N{E) — J d{E')dE' (number of levels 
up to energy E) . In contrast to our semiclassical analysis 
in the previous section, we have to take into consideration 
that our energy regime is finite, therefore we approximate 
N{E) by the generalized Weyl's law including perimeter 
and corner corrections |^. We obtain for the rectangle 
with Dirichlet boundary conditions in Fig. |^(b) 



N{E) = E 



Ij. -\- '2ly 

in 



E+± 
16 



(11) 



and for the rectangle with mixed boundary conditions in 
Fig. ic) 



NiE).E-^-^VE+^^ 



(12) 



By construction, the unfolded spectra {En} have unit 
mean level spacing. Henceforth, the tilde ~ will be sup- 
pressed. 



be observed in the magnification. The fluctuations de- 
crease with increasing energy, and they are larger than 
the statistical fluctuations < 0.5/ ^/W ~ 0.0035 due to 
the finite width W — 20 000 of the energy windows. For 
the Veech billiard, we find a slight tendency towards the 
Wigner surmise for medium energies and a slight ten- 
dency towards the Poisson distribution for high energies. 
The fluctuations in the non- Veech case are of the same 
magnitude but without clear tendency towards Wigner 
surmise or Poisson distribution. 
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A. Nearest-neighbor spacing distributions 

An important statistical quantity measuring short- 
range level correlations is the nearest-neighbor spacing 
distribution. It is defined as the probability density of 
the spacing s between adjacent levels 



1 " 

P[s) = lim - y S{s - E.,+i + E,) 



(13) 



We will compute its integral, the cumulative spacing dis- 
tribution 



P{s')ds' 



(14) 



For Poisson statistics Pp{s) = exp (— s) and /p(s) = 
1 — exp (— s), the GOE is well described by the Wigner 
surmise P\j^{s) — (7r/2)s exp (— 7rs^/4) and /vy(s) 
1 — exp (— 7rs^/4), and for the SP statistics 



ad /w ( 



Psp(s) = 4se-2-, /sp(s) = 1 - {2s + l)e 



-2s 



(15) 



Pgp(s) shows a linear increase at small s (level repulsion) 
like the Wigner surmise and an exponential fall-off at 
large s like Poisson statistics. 

In Fig. H one sees that the cumulative spacing distri- 
bution is in good agreement with the SP statistics in 
both energy regimes (the medium-energy behavior of the 
non- Veech billiard is not shown since it is similar to the 
Veech case). However, small fluctuations around SP can 



FIG. 2: Difference between the cumulative spacing distribution of 
the pure barrier-billiard levels. Below: magnification. 

The fluctuations for the Veech barrier billiard are very 
similar (but by a factor 2.5 smaller in the high-energy 
regime) than those found in right triangles [g2|. In p2[ , 
increasing fluctuations have been reported for very high 
energies above the 4 000 000th level. These fluctuations 
have been interpreted as deviations from SP leading to 
the conclusion that SP is asymptotically not the relevant 
statistics for pseudointegrable systems. In the following 
paragraphs, however, we will show that this interpreta- 
tion is doubtful. 

Let us construct an artifical SP distributed sequence 
of numbers. Take the levels of the simple rectangle in 
Fig- 0(b) given by 



(16) 



with m, n = 1, 2, 3, . . .. It has been demonstrated numer- 
ically that the nearest-neighbor spacing distribution and 
some other statistical properties of such a sequence are 
asymptotically extremely well described by the Poisson 
statistics ; see also . After ordering the levels ac- 
cording to increasing energy and removing every other 
level, the nearest-neighbor spacing distribution of the se- 
quence thus obtained obeys SP statistics [^. Figure ^ 
shows the corresponding cumulative spacing distribution 
computed numerically from 20 000 levels in three differ- 
ent regimes. The medium- and high-energy regime are 
defined as before, whereas the very-high-energy regime 
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starts at the 4 000 000th level as in Ref. |2|l. We observe 
small fluctuations around SP which decrease with in- 
creasing energy. In the medium- and high-energy regime, 
the fluctuations are of the same order of magnitude as for 
the pure barrier-billiard levels; cf. Fig. |^. We note that 
the same fluctuations are also present when I{s) ~ Ip{s) 
is plotted for the Poisson sequence given by Eq. (|l6|). 

The statistical fluctuations depend on the number of 
levels under consideration. This carries over to the to- 
tal fluctuations as illustrated for the very-high-energy 
regime m Fig. I with 10 000 and 20 000 levels, respec- 
tively. Hence, one should not compare the statistics of 
sequences with different number of levels as it has been 
done in H. 

Following the same reasoning as described above we 
have also constructed a SP sequence of 20 000 numbers 
using a conventional pseudo-random number generator. 
This reproduces the expected statistical fluctuations of 
order 0.0035. To summarize, from the fluctuations found 
numerically here and in ||2l|, ^ it is not justified to ex- 
clude SP as correct statistics for generic pseudointegrable 
systems. 
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FIG. 3: Difference between the cumulative spacing distribution of 
the artifical SP sequence. 



We mention that the distribution of spacings between 
neighboring eigenvalues of the S-matrix in an open ver- 
sion of the barrier billiard also resembles the SP re- 
sult; although the agreement is not as good as here. 



B. Next-to-nearest spacing distributions 

In the previous subsection we have seen that the 
nearest-neighbor distributions are close to the SP pre- 
diction in Eq. (]l5|). However, Eq. ( ]l5| ) is also valid for 
IP statistics. In order to distinguish between IP and SP 
statistics one has to consider other correlation functions. 
First, we choose the next-to- nearest spacing distribution 
(second-neighbor-spacing distribution) and its integral. 



For the SP statistics 
^SP(2,s) ^ 



3 

1 



/sp(2,s) = 1 - ^{As^ + 6s^ + 6s + 3)e-^' . (17) 

For IP we find analytically 

i^p(2,s) = Ae-^[l-(l + s)e-^], 



/ip(2,s) = l-|-e-'[e"^(3-|-2s)-4] 



(18) 



Figure |j shows that the cumulative next-to-nearest 
spacing distribution is in agreement with the SP statis- 
tics but not with IP statistics. Note that the fluctuations 
are two times larger as in the case of the nearest- neighbor 
spacing distribution in Fig. ^. 
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FIG. 4: Difference between the cumulative next-to-nearcst spacing 
distribution of the pure barrier-billiard levels. 

We have also investigated nth-neighbor spacing distri- 
butions P{n,s) with n = 3,4,5. Again, the distribu- 
tions differ significantly from IP statistics and are well 
described by SP statistics, even though the fluctuations 
increase slightly. A detailed discussion is left out since 
we will study long-range level correlations in a more com- 
prehensive way in the next subsection. 



Number variance 



The number variance 

nL)^{{n{L,E)-Lf) 



(19) 



is the local variance of the number n{L^ E) ~ N{E - 
L/2) — N{E — L/2) of energy levels in the interval [E 
L/2,E + L/2]. SP statistics gives §, 0, ^ 



Ssp(i) = § + ^(l-e-^^) 



(20) 



For IP statistics we get analytically a different result 

Eip(L) = L-i + (L + i)e-2^. (21) 
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Figure || reveals a substantial difference to SP for cor- 
relation lengths L > 4 in the medium-energy regime. In 
the high-energy regime the difference is smaller. Note 
that the L regime in Fig. ^ is well below the crossover 
region where the number variance begins to saturate at 
a value determined by the shortest periodic orbit ||3^ . In 
the region of large L, the number variance is related to 
the form factor (see, e.g., |p3|) by means of 



X(0) 



lim 

L — >oo 



(22) 



Using this relation we get for the Veech case K{Q) « 0.27 
at medium energies and « 0.34 at high energies (w 0.36 
for the non- Veech case). However, we do not interpret 
this result as deviation from SP since we know from Sec- 
tion in that in the Veech case if (0) does converge to the 
SP result 1/2. Hence, we conclude that the convergence 
to a stationary limit is extremely slow. The slow conver- 
gence of the spectral statistics is shared by related sys- 
tems such as right triangular billiards |2^ psf, rect- 
angular billiards with magnetic flux lines |23|, [43[, and 
parabolic maps with spin p4| . To overcome the problem 
of slow convergence, we have tried to use an extrapola- 
tion procedure described in ]2^ . However, in our case it 
does not give satisfactory results and therefore a detailed 
discussion is omitted. 
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FIG. 5: Number variance S(L) of the pure barrier-billiard levels. 



Figure |6| shows the number variance for the artifical 
SP sequence constructed from the levels of the integrable 
rectangle. The convergence in direction towards SP is 
similar, even though a bit faster, as for the pure barrier- 
billiard levels plotted in Fig. |^. In the regime of very 
high energies, the number variance is hard to distinguish 
from the SP curve. 
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FIG. 6: Number variance S(L) of the artifical SP sequence. 

D. The form factor 

The form factor K{t) can be approximated numeri- 
cally by (see, e.g., Ref. |l|) 



n 



l+n 



(23) 



We here consider only the high-energy regime, i.e. I — 
400 000 and n = 20 000. We smooth K{t; n) by averaging 
over small intervals of size Ar = 0.006. Nevertheless, 
the numerical data is quite irregular as can be seen in 
Fig. for the Veech billiard (for the non- Veech case the 
picture looks very similar). It is difficult to estimate K{0) 
directly from such kind of data, but it is clear that K{0) is 
well below the SP prediction 1/2, which is consistent with 
our former numerical results on the number variance. 

A more elegant way to compare the form factor to SP 
statistics is described in Ref. Fit K{T;n) to the 

function 



if fit (r) 



i2 - 2a + 4^V2 



(24) 



Expression (24) is the SP form factor when a = 4. There- 
fore, the quantity iffit(O) — 1/2 measures the difference 
to SP statistics. Note that i<'fit(0) in general differs from 
K{Q;n) since it depends also on K{T;n) with r > 0. 
Fitting Eq. (|2^) to our smoothed data over the range 
< T < 3, we find remarkable agreement with SP statis- 
tics: if fit (0) = 0.504 for the Veech bilhard (see Fig. |) 
and if fit (0) = 0.498 for the non- Veech bilhard. 



E. Level dynamics 

We here investigate the dependence of the energy levels 
on the system parameter l/ly. This so-called "level dy- 
namics" has been intensively studied for classically inte- 
grable and chaotic systems; see e.g. [Esl. To the author's 




FIG. 7: The form factor {g3|) of thermre barrier- billiard levels; 
I = ly/2. The smooth curve is the fit @ with a = 4.032. 

knowledge, only one pseudointegrable (for all parame- 
ter values) example, the "square torus billiard" ^ and 
a generalized version of it |}45| , has been studied in this 
regard. 

A typical situation for the pure barrier-billiard levels 
is displayed in Fig. ^. The global increase of the levels 
(not unfolded) is due to the fact that the smooth part of 
the integrated density of states in Eq. (|l2|) decreases as 
I is increased. Apart from this rather trivial fact we ob- 
serve a number of interesting features: (i) the levels tend 
to avoid each other. Closer examination of the available 
numerical data indicates that there are no level crossings. 
That means for fixed parameter value there are no degen- 
eracies in the spectrum, which is consistent with the SP 
and the GOE prediction for the nearest-neighbor spac- 
ing distribution P(0) = in agreement with our former 
numerical results. The total absence of level crossings is 
in contrast to the situation in the square torus billiard 
where crossings can appear for parameter values at which 
the billiard is almost integrable 0]. (ii) Large areas free 
of levels exist similar as in integrable systems and dif- 
ferent to fully chaotic systems. This is consistent with 
Poisson and SP statistics which both predict a slower 
fall-off of P{s) at large s than GOE statistics does, (iii) 
There exists an unusual structure of plateaus interrupted 
by steep segments not only near avoided crossings but 
also fairly far away from avoided crossings. Observation 
of the energy eigenfunctions reveals that plateaus (steep 
segments) correspond to parameter values at which the 
corresponding eigenfunction has small (large) amplitude 
at the upper end of the barrier. Hence, varying the bar- 
rier length has no (strong) influence on the wave pattern 
and on the energy, resulting in a plateau (steep segment). 

Interestingly, the abrupt changes in the slopes fairly 
far away from avoided crossings can be simulated in a 
natural way by an artifical SP sequence constructed as 
described before by removing every other level of an in- 
tegrable system. Two neighboring levels of an integrable 
system typically cross each other when a parameter is 
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FIG. 8: Pure barrier-billiard levels as functions of the barrier 
length Omiy < I < 0.95ly. 



varied as sketched in Fig. Q Removing the second level 
(measured from below for each value of the parameter) 
gives the solid, nondifferentiable line. This could produce 
the kind of abrupt changes seen in Fig. ^. Of course, 
the slope of finite-energy levels cannot change discontin- 
uously. Real discontinuouities can only be expected in 
the semiclassical limit. 



E 




parameter 



FIG. 9: Sketch of the local level dynamics of the artifical SP se- 
quence. 



V. CONCLUSION 

In this paper we have studied the energy levels of 
pseudointegrable barrier billiards. Focusing on the pure 
barrier-billiard states, we have found numerically that 
the nearest-neighbor spacing distributions and next-to- 
nearest spacing distributions agree with the semi-Poisson 
(SP) statistics which is obtained by dropping every other 
number from a random sequence. The number variance 
and the spectral form factor agree with SP, even though 
long-range correlations seem to converge rather slowly. 
Moreover, the level dynamics is consistent with SP statis- 
tics. Even though we have considered an high-energy 
window (20 000 levels starting at the 400 000th level) we 
cannot exclude that at larger energies a different scenario 
takes place. However, our analytical result for the spec- 
tral form factor for a Veech barrier billiard, K{t) ^1/2 
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as T — > 0, gives us some confidence that the spectral 
statistics of barrier biUiards are indeed close to SP. 

Due to the slow convergence of the spectral statistics in 
polygonal billiards, and other diffractive systems, semi- 
classical methods as shown here and in |23 have to be 
extended in the future to higher order in r (as in ||4^ 
for rectangular billiards with point-like singularities) and 
to other polygons in order to clarify the role of the SP 
statistics in pseudointegrable systems. 
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